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ABSTRACT 

Little attention has so far been paid to the division of the observed population between 
pulsars of the two spin directions that are possible. Almost all pulsars with positive 
corotational charge density at the polar caps are expected to satisfy space-charge lim- 
ited flow boundary conditions. Charge separation by blackbody photo-electric tran- 
sitions in moving ions limits the acceleration potential, analogously with the more 
usually considered pair creation. But the limitation is more severe so that proton and 
ion energies can be relativistic but not ultra-relativistic, and these allow the growth of 
Langmuir-mode induced turbulence that couples directly with the radiation field, as 
shown by Asseo, Pelletier & Sol. The consequences of this, and of the several possible 
physical states of the polar cap, are described, qualitatively, as possible explanations 
for the complex phenomena of nulls, subpulse drift and mode-switching observed in 
sub-sets of pulsars. 
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1 INTRODUCTION 

The past twenty years have seen a number of developments 
that are likely to be relevant to pulsar radio and X-ray emis- 
sion. Apart from the many improvements in the quality and 
scope of observations, there has been the recognition of the 
• • role of the Lense-Thirring effect in acceleration at the polar 
. ' ca P (Muslimov & Tsygan 1992) and of the significance of 
k> plasma turbulence (see the review of Melrose 2000). Most of 
< I , the very many theoretical papers published have assumed a 
C3 ■ relative orientation of rotation spin fl and polar-cap mag- 
netic flux density B such that tt ■ B > with Goldreich- 
Julian charge density paj < and electron acceleration. In 
this case, there can be no doubt that the electron work func- 
tion is so small that the space-charge limited flow boundary 
condition on the corot at ing- frame electric field, Em = 0, is 
satisfied on the polar-cap surface at all instants of time. 
(The subscripts || and _L refer to directions locally parallel 
with and perpendicular to B.) There appears to be no rea- 
son why neutron stars with ft ■ B < and paj > should 
not be present in the pulsar population, but in this case, 
the best existing calculations of ion separation energies (see 
Medin & Lai 2006) combined with estimates of the polar-cap 
temperatures can be used to obtain values of the minimum 
magnetic field for which the boundary condition Eii ^ is 
possible. These are typically of the order of 10 14 G (Jones 
2011) and are indirect evidence that all except possibly a 
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very small number of pulsars with spin direction such that 
tt ■ B < satisfy the space-charge limited flow boundary 
condition rather than the Ey 7^ condition assumed in the 
classic polar-cap model of Ruderman & Sutherland (1975). 

The concept of localized regions of electron-positron 
pair production, referred to as sparks, moving over the 
neutron-star surface within the polar cap in an organized 
way, has been adapted from the Ruderman-Sutherland 
model and the Em 7^ boundary condition, but is now 
widely used phenomenologically, with no reference to bound- 
ary conditions or to the sign of ft ■ B, in describing observa- 
tions such as nulls and subpulse drift. In the tt ■ B > case, 
with space-charge limited flow, no physical basis for such a 
phenomenological description has so far been published, the 
reason being that electrons are the only possible negatively- 
charged particles for outward acceleration and the space- 
charge limited flow boundary condition Ey = can be satis- 
fied at any conceivable polar-cap surface temperature. The 
state of the accelerated plasma is then simply a problem in 
electromagnetism with defined boundary conditions. How- 
ever, for CI ■ B < 0, it has been shown that subpulse forma- 
tion and movement, either organized or disorganized, over 
the polar-cap surface is quite natural and should be observ- 
able in most radio pulsars with this spin direction (Jones 
2010a, 2011; hereafter Papers I and II). The basis for this 
conclusion is an examination of the physical processes oc- 
curring in electromagnetic showers formed at the polar-cap 
surface by photo-electrons accelerated inward. The problem 
of finding the state of the accelerated plasma is not restricted 
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to electromagnetism but also depends on proton production 
by the formation and decay of the giant dipole resonance in 
the electromagnetic showers. 

Subpulse modulation (the wide variations of intensity 
at a fixed longitude in a sequence of observed pulses) has 
been recognized as an almost universal property of pulsars 
ever since their discovery. Weltevrede, Edwards & Stappers 
(2006) found it present in 170 out of a sample of 187 pulsars 
selected only by signal-to-noise ratio, evidence that its pres- 
ence is very likely to be independent of the sign of CI ■ B. 
Strong plasma turbulence is now viewed as the most proba- 
ble source of radio emission (see Melrose 2000). Its develop- 
ment from a quasi-longitudinal Langmuir mode and the pos- 
sible formation of a random array of stable Langmuir soli- 
tons has been investigated by Asseo, Pelletier & Sol 1990, 
and by Asseo & Porzio 2006. It has been shown recently 
(Jones 2012; hereafter Paper III) that the formation of this 
mode is also possible in a two-component beam consisting 
of protons and relativistic but not ultra-relativistic ions that 
can be formed as a consequence of the reverse flow of photo- 
electrons to the polar cap. Thus it is not surprising that 
subpulse modulation is almost universal. That a secondary 
low-energy electron-positron plasma is not the only possi- 
ble source of plasma turbulence is also consistent with the 
fact that common characteristics of coherent radio emission 
are observed in the pulsar population even though the in- 
ferred polar-cap magnetic flux densities vary by five orders 
of magnitude. Pulsars with either sign of CI ■ B are there- 
fore expected to produce individual pulses of similar form 
although if there is no pair production, there may be differ- 
ences that are difficult to detect such as in the distributions 
of spectral index caused by the different plasma frequencies 
and Lorentz factors in the two cases. But the present paper 
(see also Paper II) proposes that phenomena such as nulls, 
mode-changes and organized sub-pulse motion are restricted 
to the CI ■ B < case. These involve time-scales that have 
no obvious explanation in the CI ■ B > direction of spin. 

Existing observations show that subpulse drift is a fre- 
quent though not universal phenomenon. Weltevrede, Ed- 
wards & Stappers found that 68 of their sample of 187 pul- 
sars show observable subpulse drift. There appears to be no 
published survey either of nulls or of mode-changes based 
on a carefully defined sample as in the above work. But the 
most recent and largest tabulation of data (Wang, Manch- 
ester & Johnston 2007) shows that nulls must be a moder- 
ately frequent phenomenon. Systematic information about 
mode-changing is much more sparse (see Kramer et al 2006) . 

Models of subpulse formation and drift recently devel- 
oped (see, for example, Gil & Sendyk 2000) are based on the 
classic paper of Ruderman & Sutherland (1975) and assume 
the f2 • B < case with the polar-cap surface boundary 
condition Ey 7^ 0. These have all the attributes of a good 
physical model, particularly simplicity and predictive capa- 
bility in that it is possible to write down an expression for 
the E x B subpulse drift velocity. But the boundary condi- 
tion requires surface magnetic flux densities approximately 
two orders of magnitude larger than the dipole fields inferred 
from the spin-down rates. Although the presence of highly 
complex field geometries formed from higher-multipole com- 
ponents cannot be excluded, it would be strange if they were 
present with the frequency indicated by the survey of Wel- 
tevrede, Edwards & Stappers. Consequently, Papers I and II 



investigated physical processes at the polar-cap surface and 
have shown that the composition of the plasma accelerated 
under the space-charge limited flow (SCLF) boundary condi- 
tion is usually neither time-independent nor uniform over its 
whole area. Localized areas in which the accelerated plasma 
is suitable for the growth of the quasi-longitudinal Langmuir 
mode exist naturally and can move in either an organized 
or disorganized way, so acting as sources for subpulses with 
the same properties. Thus plasma systems analogous with 
Ruderman-Sutherland sparks are present under the En = 
boundary condition. But we emphasize that their motion is 
unrelated with E x B drift but is determined by the time 
t p in which electromagnetic shower protons diffuse to the 
surface. 

This paper attempts to present a description of the early 
stages of the coherent emission mechanism that is very dif- 
ferent from the canonical. In particular, the unexpected con- 
clusion of Paper III that relatively low-energy ion and proton 
beams can be formed means that it is necessary to re-assess 
the model of sub-pulse drift discussed in Paper II. Polar caps 
of pulsars with spin CI ■ B < are seen to be very complex 
systems compared with those of the electron acceleration 
CI ■ B > case. But we must presume that both are present 
in the neutron-star population and so it is worth examining 
how the properties of the observed pulsar population may 
be related with spin direction. 

Before proceeding further, it must be admitted that the 
work has at least one drawback. Processes of the kind consid- 
ered in Papers I and II, occurring at a continuously chang- 
ing real condensed-matter surface, are not necessarily sus- 
ceptible to the formulation of simple physical theories. The 
behaviour of such systems can be varied and complex. But 
this is not inconsistent with the observed properties of many 
individual pulsars. 

Section 2 summarizes the assumptions about the prop- 
erties of the open region of the magnetosphere that we are 
obliged to make in both qualitative and quantitative work. 
Strong plasma turbulence is possible only for limited inter- 
vals of particle energy at altitudes 2 above the polar-cap 
surface smaller than z ~ 107?, where R is the neutron-star 
radius. This is described in more detail in Section 3. Sec- 
tion 4 draws on the work of Paper III to obtain values of 
the parameter K , which is the number of protons produced 
by photo-electrons per unit nuclear charge accelerated. We 
then give in Section 5 a description of the polar cap showing 
how subpulses are formed and that either chaotic or orga- 
nized motion is possible. An attempt is made to show that 
this model of the polar cap provides a physical basis for 
phenomena such as nulls, subpulse drift, and mode changes. 



2 THE OPEN MAGNETOSPHERE AT THE 
POLAR CAP 

At the base of the magnetosphere there is an atmosphere 
of ions, assumed to be in local thermodynamic equilibrium 
(LTE), having a scale height of the order of 10 _1 cm. Its 
total mass is temperature-dependent but is broadly equiva- 
lent to 10 _1 — 10 1 radiation lengths. Thus it contains only 
the earlier part of the electromagnetic showers formed by 
the reverse-electron flux. There is fractionation of the ionic 
charge to mass ratio. Ions with the highest value are con- 
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centrated at the top of the atmosphere and so preferentially 
enter the acceleration region. The proton number density 
is many orders of magnitude smaller than that of the ions 
and consequently has no effect on the structure of the atmo- 
sphere. Protons are not in static equilibrium within the LTE 
ion atmosphere but are subject to an outward-directed net 
gravitational and electrostatic force. Their motion is there- 
fore a combination of diffusion and drift with a characteristic 
time t p elapsing between formation and arrival at the top 
of the atmosphere. If the rate of formation is such that the 
Goldreich- Julian flux would be exceeded, the excess protons 
form a thin atmosphere in equilibrium at the top of the 
LTE ion atmosphere. We refer to Sections 2, 3 and 5 of Pa- 
per II for a more complete account of these topics. Values 
r p ~ 1CP 1 — 10° s are estimated. 

The polar-cap radius adopted here and in Paper III is 
based on dipole-field geometry and is that given by Harding 
& Muslimov (2001), 

for rotation period P, in which /(l) = 1.368 for the neutron 
star of mass IAMq and radius R = 1.2 x 10 6 cm. As in Paper 
III, we use the fact that the open magnetic flux lines are 
contained within a narrow tube, assumed to be of circular 
cross-section whose radius, uo(z) = iio(0)77 3,/2 for a dipole 
field, increases only very slowly as a function of altitude on 
the magnetic axis z = (rj — 1)R at radii r\ ~ 10 or smaller. 
Therefore, we can adequately represent a small section of the 
tube at altitude z as a right cylinder of locally-defined radius 
u = Uo(z). Then the electrostatic potential $ at z > no(0) 
and under the SCLF boundary conditions is approximately 



$(«, Z) = 7T («o(z) - M 2 ) (p(z) - PG.j(z)) , 



(2) 



in terms of the difference between the charge density p and 
the Goldreich- Julian charge density pcj which is almost con- 
stant over the cross-sectional area of the tube. We refer again 
to Harding & Muslimov (2001) for these quantities at radial 
coordinate 77, 



PGJ 



and, 



P = 



Bf(y) 
cPar) 3 f(l) 



COS 1p 



Bf(v) 



(1 — k) cos ip, 



(3) 



(4) 



cParj 3 f(l) 

in which B is the surface magnetic flux density. Equation (4) 
gives the SCLF charge density undisturbed by photo-electric 
transitions or electron-positron pair creation. In these ex- 
pressions, k is the dimensionless Lense-Thirring factor which 
we assume to be k — 0.15, and tp is the angle between fi 
and B. Equation (4) also assumes that particle velocities, in- 
cluding those of protons and ions, differ negligibly from the 
velocity of light. The function f(r])/af(l) including the red- 
shift factor a is a slowly- varying function of 77 and is approxi- 
mately unity at low altitudes. Thus pcj(0) = (1 — k)pgj(oo) 
where poj(oo) is the flat space-time Goldreich- Julian charge 
density. Equations (3) and (4) exclude terms that can be- 
come significant at higher altitudes for which our geometri- 
cal model also fails, but nevertheless equation (2) is adequate 
for the estimation of proton and ion energies and of photo- 
electric transition rates. It is convenient to use the energy 



units of particle physics for polar-cap processes and there- 
fore we define a potential energy difference V(u, z) = — e$. 
It has the value V ma x (u, z) on a flux line with polar-cap 
coordinate u for the undisturbed SCLF charge density of 
equation (4). It is also useful to have a simple expression for 
Vmax on the magnetic axis, 



V m ax(0, CO) 



2n 2 R A KeB 
c 2 /(l)P 2 



= 1.25 x 10 d 



1 Bn 



GeV, 



(5) 



for acceleration above the polar cap, where B12 is the surface 
magnetic flux density in units of 10 G. 

The above statement is subject to a serious qualifica- 
tion concerning the definition of the open magnetosphere. 
If many secondary low-energy electron-positron pairs were 
created per primary particle accelerated, there would be no 
difficulty in assuming the unconstrained outward flow of pri- 
mary particles and secondary pairs on any magnetic flux line 
intersecting the light cylinder. The separate velocity distri- 
butions of secondary electrons and positrons would have the 
capacity to adjust so as to maintain the condition Em = 
at all points in the open region. But there must be serious 
doubts about the universal existence, in the pulsar popu- 
lation, of the necessary density of secondary-pair plasma. 
Apart from fl ■ B < pulsars in which there may be no pair 
creation, calculations of the secondary pair density produced 
by the inverse Compton scattering (ICS) of polar-cap pho- 
tons by outward accelerated electrons in the ft ■ B > case 
(Hibschman & Arons 2001; Harding & Muslimov 2002) pre- 
dict in many pulsars pair densities that are small compared 
with unity, per primary electron. It is true that more com- 
plex field geometries, such as the off-set dipole, are possible 
and may enhance pair production (see Harding & Muslimov 
2011) but there remain at least some hints of a problem in 
the pair density. In the absence of a pair plasma with the re- 
quired properties, motion in the magnetosphere beyond the 
77 ~ 10 region depends on the sense of curvature of the open 
flux line: there is further acceleration for | cos^>| increasing, 
but deceleration in the opposite case which is capable of can- 
celling the acceleration at lower altitudes. A section of the 
polar cap is then a dead zone from which there is no outward 
particle flux. Whilst the change in active polar-cap shape is 
not very significant, the potential given by equation (2) is 
reduced, very roughly, in proportion to the decrease in area. 
The component Ey in the dead volume is very small and 
is that required to maintain a charge-separated equilibrium 
with gravity and the centrifugal force. 

A further source of uncertainty is the consistency with 
which the total electric charge of the neutron star remains 
constant. The rate of loss of charge at the polar caps is 
q — 27ruo(0) pgj (0)c whose compensation is likely to be easy 
if there is strong secondary pair creation, but is otherwise ob- 
scure. An order of magnitude estimate of the relevant time- 
scale can be found by dividing a typical light-cylinder field 
strength by the rate of change of the light-cylinder radial 
electric field given by q. This is, 



BR 3 
Rlc 



2P 



(6) 



where Rlc is the light-cylinder radius. This indicates that 
if there is to be no substantial change in the structure of 
the magnetosphere, cessation of the processes maintaining 
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charge constancy is possible only for time intervals no more 
than of the order of P. 

But the great uncertainty affecting fl ■ B < pulsars 
is in the atomic number of matter at the neutron-star sur- 
face. If this were small (Z ~ 6) there would be a negligible 
reverse-electron flux and, as we shall see in Section 3, no 
coherent radio emission unless curvature radiation electron- 
positron pair production is possible. There is some indirect 
observational evidence for Z — 6 (Ho & Heinke 2009; see 
also Bogdanov & Grindlay 2009, Yakovlev et al 2011, and 
Zhu et al 2009) but it does not appear strong. Theoret- 
ical work does indicate the possibility of atomic numbers 
other than the canonical Z = 26. Calculations by Hoffman 
& Heyl (2009) for isolated neutron stars with no mass ac- 
cretion show the presence of Z — 14, though with small 
total mass, but Pearson, Goriely & Chamel (2011) find a 
range of values Z > 26. The mass per unit area removed 
by a Goldreich- Julian current density of nuclei with mass 
number A in an fl ■ B < pulsar in a life-time t is, 

Am p f —Bcostjj , — 2Am p Bt cost/) . , 

— J -7Pl¥T dt * zTFW) ' (7) 

assuming that B and ip are time-independent. This is ~ 
4 x 10 n Bi2P _1 g cm" 2 at t = 10 6 yr, and over the whole 
surface of the star amounts to a layer of no more than 4 x 
10 _9 -Bi2-P _1 Mq. It is a quantity that could easily be the 
result of unknown details of the formation process such as 
fall-back. In view of this, direct experimental evidence of the 
nuclei that are actually present at the polar-cap surface is 
clearly needed. 

The uncertainty in the atomic number of polar-cap sur- 
face nuclei has an obvious effect on photo-electric transi- 
tions, and therefore on the reverse-electron flux and the ac- 
celeration potential difference for particles moving on any 
given magnetic flux line. The dependence of 3> on p(z) 
in equation (2) means that the production of oppositely 
charged particles in a charge-separated region at altitude 
z under the SCLF boundary condition tends to cancel the 
local electric field Eii . This is well-known in the case of 
electron-positron pairs but is, of course, also true for photo- 
electrons (see Paper III). 

Equation (2) is not valid at low altitudes z ~ tto(0). 
Near the polar-cap surface, z <C «o(0), Em is the inertial ac- 
celeration field first described by Michel (1974). The problem 
of finding the electric field is one-dimensional in z and the 
inertial Ey is not large. Its contribution to particle acceler- 
ation is small compared with that of the Muslimov-Tsygan 
effect. We refer to Paper III for an approximate model for 
this region. 

The polar-cap and whole-surface neutron-star temper- 
atures are treated as distinct in Paper III and in the present 
work. Calculated photo-electric transition rates for polar- 
cap blackbody photons are negligibly small above an alti- 
tude z = h « 4tto(0), principally because the photon mo- 
mentum component perpendicular to B becomes small at 
higher altitudes. This is unfavourable both for the Lorentz 
transformation to the rest-frame of the accelerated ion and 
for the cross-section to the lowest Landau state of the emit- 
ted electron (see Paper III). 

The SCLF boundary conditions lead to an electric field 
having a complex form in the z < h region which is un- 
likely to be well represented by equation (2) . There also exist 



the complications, just discussed, of the inertial En and of 
photo-electric transitions in this interval. But the boundary 
conditions also limit the extent to which p(z) can change 
if the condition En > is to be satisfied. Thus the differ- 
ence between the mean ion charges Zh at z — h and Z in 
the assumed local thermodynamic equilibrium of the atmo- 
sphere at z — cannot be large. Adjustments in p(z) and in 
Eii maintain this state by limiting ion Lorentz factors and 
hence photo-electric transition rates. The ratio, 

p(h) = KZ S + 2Z h - Z x 
P(0) ~ KZ S +2Z-Z 00 1 

including ion, proton and electron charge densities is always 
very nearly unity at z < h. Here K is the mean number of 
protons produced in electromagnetic showers per unit nu- 
clear charge accelerated. The mean charge of nuclei reaching 
the surface after moving through the electromagnetic shower 
region is Z s , and Zoo is the ion charge at z max , where photo- 
electric transition rates become negligibly small. Because the 
acceleration in this interval is minor compared with that at 
z > h we can assume that the effect of this region is to in- 
ject ions of mean charge Zh and Lorentz factor ju into the 
main acceleration region z > h in which p(z) is modified by 
photo-electric transitions and possibly by electron-positron 
pair creation. Although the mean value of Zh — Z is small, 
it is undoubtedly finite, and the associated reverse-electron 
flux gives a modest energy input, , to the polar-cap sur- 
face even in the absence of any whole-surface contribution 
to photo-ionization at z > h. We refer to Paper III for its 
approximate size, ~ 20 GeV. 

For a relatively small fraction of pulsars, pair creation 
by the conversion of curvature radiation (CR) photons is 
possible. But for most fl • B > pulsars, it is believed that 
the source of pairs is the conversion of outward-directed ICS 
polar-cap photons (see Hibschman & Arons 2001, Harding 
& Muslimov 2002). However, in the case of the fl ■ B < 
pulsars considered here, there is an essential difference in 
that the electron-photon momentum in the frame of the ro- 
tating neutron star is directed inwards. There is likely to be 
some pair creation at low values of z from 7-rays emitted by 
the capture of shower-generated neutrons at the polar-cap 
surface, but the rate is difficult to estimate. Thus pair cre- 
ation rates in most fl ■ B < pulsars are likely to be one 
or two orders of magnitude smaller than for those with spin 
direction fl ■ B > 0. But even for the latter spin direction, 
the number of ICS pairs formed per primary electron can 
be much smaller than unity (see Harding & Muslimov 2002, 
also Fig. 8 of Hibschman & Arons). Thus for the fl ■ B < 
spin direction, there is a very real problem in understand- 
ing how a density of low-energy secondary electron-positron 
pairs adequate for coherent radio emission can be formed. 



3 CONDITIONS FOR STRONG TURBULENCE 

There is now a consensus that, in almost all observable ra- 
dio pulsars, the source of the coherent emission lies within 
radii r\ ~ 10 1 — 10 2 and its energy is derived from parti- 
cle acceleration at the polar cap. There is also a growing 
consensus that the formation of strong plasma turbulence is 
involved (Melrose 2000). The unstable mode that develops 
into strong plasma turbulence and transfers energy from par- 
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tide beams to the electromagnetic field must therefore have 
a high growth rate. This severely constrains the energy of 
the beam particles. The equations of motion for the system 
are the Maxwell equations and those for a relativistic parti- 
cle fluid. In a strong magnetic field, the mass-to-charge ratio 
present in the equations is niijfq' 1 for particles of species 
i, where qi is the charge and 7, the Lorentz factor. The lon- 
gitudinal effective mass is nnjf and this factor necessarily 
determines the growth rate of any unstable mode. 

We consider in this paper, specifically, the growth of 
a quasi-longitudinal Langmuir mode (Asseo, Pelletier & Sol 
1990) from which a random array of stable Langmuir solitons 
may develop (Asseo & Porzio 2006). A different longitudinal 
mode has been studied by Weatherall (1997, 1998) and may 
be the source of the nanosecond pulses observed in certain 
pulsars (Hankins et al 2003; Soglasnov et al 2004). (We re- 
fer to the paper of Asseo & Porzio for an account of the 
history of plasma turbulence in relation to pulsar physics.) 
The Langmuir mode and solitons have components Em 7^ 0, 
Ex 7^ and Bt perpendicular to Ex and so can transfer 
energy to the radiation field directly (see Asseo, Pelletier & 
Sol 1990). Hence the direct transfer of energy to the radi- 
ation field is determined only by the longitudinal effective 
mass and by the departure of the mode wave-vector from 
the purely longitudinal state. Examination of growth rates 
shows at once that primary electrons or positrons of ~ 10 3 
GeV energy have longitudinal effective masses so large that 
they are unable to participate in the mode. The beams that 
interact in fl ■ B > pulsars can be only secondary electrons 
and positrons whose velocity distributions satisfy the rela- 
tivistic Penrose condition (Buschauer & Benford 1977). But 
ion and proton beams that have longitudinal effective masses 
of the same order as secondary electrons and positrons can 
be produced in fl ■ B < neutron stars. Thus it is antici- 
pated that pulsars of both directions exist and produce very 
similar coherent emission. 

Beams of cold ions and protons (or, hypothetically, 
cold electrons and positrons) permit growth of the quasi- 
longitudinal Langmuir mode for which the dispersion rela- 
tion is, 



kl 



1 - 



(8) 



for angular frequency uj, wave-vector k, and ion and proton 
(or, electron and positron) velocities /3i and /3a, respectively. 
The remaining quantities are uj* = j^ 3 OJi in which, 

Wi = , (9) 



where qi, riii and Ni are respectively the charge, mass and 
neutron-star frame number density of each component i = 
1, 2. The transverse field components are, 

k,E« 



E 1 



B t = 



2{k\\ - 
k±E» 



k±k 2 E ll 
w(fcj? - U) 2 ) ' 



(10) 



The plasma frequency in the rest frame of a component is 
"f~ 1 ^ 2 uii and so, following Asseo et al, a natural wavenumber 
for which to seek a solution is defined in terms of the slower- 
moving component, and is kn = 2o;*7i. With the definition 



of a new variable s such that, 
u> — k\\Pl = 0J\ (1 + a), 
and, 

fen 



27? 



valid for 71 ^> 1 equation (8) can be expressed as, 
1 



1 C _ kl-fi 



(s + 1) 2 (S + M ) 2 



kh 



(11) 

(12) 
(13) 



in which fi = 71/72 and C = u>2 2 /lu* 2 . In the limit k± — t 0, 
equation (13) is a quartic with two real and two complex 
roots, whose properties were investigated in Paper III. The 
amplitude growth rate was found to be of the order of Imaj ~ 
0.2a>i and is a slowly varying function of the ratio 71 /72 and, 
for ions, the proton-ion number density ratio. But we have 
to consider the quasi-longitudinal case of k± 7^ in which 
the transverse field components are non-zero. Equation (13) 
is then less transparent, being a quintic with three real and 
two complex roots, but we can use perturbation theory to 
estimate the extent to which a non-zero k± changes a given 
complex root. Let s = so + 8s, where so is a root of the 
quartic. Then, 



Ss 



C 



(so + l) 3 (so + m) 3 



; 2 2 
kffSo 



(14) 



From equations (10), we see that \E±\ = \E\\ for k± = 
fcn |s|/ti. Then for typical values of C, \i and so, we have, 



\Ss\ 



1 2 2 
fcx7i 

kfi I s | 



|£0| 

2 ' 

7! 



(15) 



indicating that for a significant range of k± , the growth rate 
does not differ much from the k± = value. 

The growth rate decreases with distance above the polar 
cap and so it is appropriate to obtain the amplitude growth 
factor expA at a given radius r = r/R by integration. It is, 

1/2 



A = Ims 



-WirctiqiR Bcosip 
Pmic 3 7 3 



n 



-1/2 



)■ 



(16) 



assuming a dipole field and constant 71. Here, the ion (or 
electron) component number density in the frame of the ro- 
tating neutron star is a fraction ai of the Goldreich- Julian 
number density. For Ims = 0.2 and costp — —1, the factor 
is, 

1/2 



A = 2.4 x 10 3 



BiaQimai 
Pemi7 3 



(17) 



where m is the electron rest mass. We can see that the 
growth factor increases most rapidly at small rj and is sub- 
stantially developed at 77 = 2. The only particle property 
that it depends on is the charge to longitudinal effective 
mass ratio. A minimum value necessary for the development 
of nonlinearity is presumably given by the size of particle 
number fluctuations within a cubic-wavelength volume at 
the Goldreich- Julian number density. This indicates a mini- 
mum A ~ 20. We assume A = 30. For a typical pulsar with 

PB^2 = 1, equation (17) then requires 71 < I76a\^ 3 for 
1/3 

electrons or 71 < llQj for ion beams at r) — 2 increasing 
to 71 < 311«i' /3 and 71 < 20a} //3 , respectively, at 77 = 10. 
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Development of turbulence or of stable solitons is then pos- 
sible for either electrons or ions although for electrons, the 
cold beam hypothesis is clearly unrealistic and the distri- 
butions of electron and positron velocities must satisfy the 
Penrose condition. 

Quasi-longitudinal modes transfer energy directly to 
the radiation field at angular frequencies of the order of 

1/2 

ujq — 27/ wi. The actual maximum in the intensity dis- 
tribution presumably depends on the extent to which de- 
veloped turbulence moves energy toward higher wavenum- 
bers. This is unknown except that some examples of soliton 
formation considered by Asseo & Porzio have maxima at 
uj ~ 1.5wo- As a function of 77 and for cos?/; = — 1, the fre- 
quencies vq — ujo/2it are, 




respectively, for electrons and for ions of mass number A and 
charge Zoo- These values differ by two orders of magnitude 
and suggest the existence of two classes of pulsar with very 
different radiation spectra though with individual pulses of 
similar form. There have been relatively few measurements 
of flux density below 400 MHz but those of Deshpande & 
Radhakrishnan (1992) and of Malofeev et al (1994) both 
confirm that low-frequency spectral cut-offs are usually be- 
low 100 MHz and in some cases, below 50 MHz. 

There have been many calculations of pair formation 
densities above polar caps, but of particular interest is the 
paper of Harding & Muslimov (2011) who have studied the 
effect of a specific type of field, that of an off-set dipole. Pair 
formation densities, assuming outward moving ICS photons 
to be the source, are rapidly increasing functions of the de- 
gree of off-set. The most extreme case in which pair forma- 
tion might be thought doubtful is that of PSR J2 144-3933 
which has a period P — 8.51 s and whose observers (Young, 
Manchester & Johnston 1999) commented specifically on the 
problem. The value of its parameter BP~ 2 = 2.87 x 10 10 G 
s~ 2 is almost an order of magnitude smaller than the ap- 
parent cut-off value of 2.2 x 10 11 G s -2 in the distribution 
of that quantity obtained from the ATNF pulsar catalogue 
(Manchester et al 2005) representing a death-line in the P-P 
distribution below which pair creation is not possible. It is 
old (2.7 x 10 8 yr) with a relatively low mean flux density at 
400 MHz, and may be observable only as a consequence of 
close proximity (170 pc). Although the maximum accelera- 
tion potential given by equation (5) is very small, Vmax = 36 
GeV, ICS pair formation is possible in principle although a 
very high degree of off-set would be required (see Figs. 7-9 
of Harding & Muslimov 2011). However, given a spin direc- 
tion f2 ■ B < 0, accelerated ions would have Lorentz factors 
appropriate for growth of the quasi-longitudinal mode for 
any field geometry. But the value of vo given by equations 
(18) is rather small in this case, even assuming that the ra- 
diation field is largely decoupled from the magnetosphere at 
quite low altitudes 77 ~ 2. (For this spin direction, the pair 
multiplicities produced by ICS photons would probably be 
one or two orders of magnitude smaller than for J7 • B > 
because the momentum of the electron-photon system is in- 
ward directed.) But more extensive observations at low fre- 



quencies may reveal whether or not J2 144-3933 is unique or 
merely a nearby example of an otherwise unobserved popu- 
lation. 



4 INSTABILITY AND PROTON 
PRODUCTION 

Owing to our late recognition that photo-electric transitions 
have a very significant effect on particle acceleration under 
SCLF boundary conditions, it is necessary to revise the es- 
timates of the reverse-electron energy flux that were used 
in Papers I and II. Photo-electric transitions and the conse- 
quent reverse-electron flux produce changes in mean nuclear 
and ion charge that we can summarize here as follows. The 
initial nuclear charge is Z (canonical value Z — 26) but for- 
mation and decay of the giant dipole state in electromagnetic 
showers reduces it to Z s at the top of the neutron-star atmo- 
sphere. We assume local thermodynamic equilibrium (LTE) 
in this region with ion charge Z. Then ions with charge Zu 
are injected into the main acceleration region at z = h. Typ- 
ically, Zh ~ Z -\- 1 owing to interaction with polar-cap black- 
body photons. Then further acceleration produces more ex- 
tensive ionization to Zoo ^ Z s by the whole-surface black- 
body field with Z(z) as the altitude-dependent ion charge 
during this process. 

In papers I-III we defined KZ S to be the number of 
protons produced per ion of nuclear charge Z s accelerated. 
The initial problem treated here is of acceleration given by 
equations (2) - (4) with p(z) independent of time and of 
position u. Then, 

p(z) _ KZ S + 2Z(z) - Zoo 
p(h) ~ KZ S + 2Z h - Zoo ' 

and the total photo-ionization at altitude z is related to the 
potential at that point by, 

2Z(z) - 2Z h _ k(^ 3 - rT 3 ) ( V(u, z) \ 
KZ S +2Z h - Zoo 1 - ^ V m ax(u, z) ) ' 1 ; 

where r/u = 1 + h/R. This result takes into account pro- 
tons, ions and the reverse flux of photo-electrons, but not 
CR electron-positron pair creation or the conversion of ICS 
photons. 

Photo-electric transition rates are rapidly increasing 
functions of the ion Lorentz factor 7. The reverse-electron 
flux they produce leads to a charge density that increases 
with altitude, so reducing the potential given by equation 
(2). Thus we might expect a self-consistent set of functions 
Viz, u), p(z, u) and 7(2, u) that are time-independent. In or- 
der to investigate this possibility, it is necessary to calculate 
the reverse-electron energy flux taking into account the fact 
that photo-electric transitions reduce the potential energy 
V — — derived from equation (2) to levels below Vmax- 
Broadly, we expect that 7(2:) will increase from its initial 
value jh but that V(u, 00) < V m ax{u, 00), and our initial at- 
tempt to model this is to assume that V(u, z) = V m ax(u, z) 
until a cut-off V c is reached. 

With the procedures and assumptions of Paper III, we 
calculate the mean reverse-electron energy e s per ion re- 
sulting from interaction with whole-surface blackbody pho- 
tons, and the mean ion charge Zoo- Transition rates are ob- 
tained to an altitude z max = 3R, beyond which they become 
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small and equation (2) predicts negligible further accelera- 
tion. The polar-cap temperature is treated as distinct from 
the general surface temperature of the star and its value is 
equal to that generated by the reverse-electron flux which 
produces protons at the Goldreich- Julian current density. 
A flux exceeding this produces an accumulation of protons 
at the top of the neutron-star atmosphere which, until ex- 
hausted, reduces the ion component of the current density 
to zero. We adopt equation (33) of Paper I, with cos^ = — 1 
and rotation period P = 1 s, and a typical proton produc- 
tion rate per unit shower energy of W p mc 2 — 2 x 10 -4 , that 
is, one proton per 2.5 GeV shower energy. This quantity is a 
slowly varying function of B and Z and has been estimated 
in Papers I and II. Then (eh + e s )W p = KZ B with eh ~ 20 
GeV being the contribution of polar-cap blackbody photons. 
The temperatures are T pc = 0.74 x 10 6 K, 0.82 x 10 6 K and 
1.10 x 10 6 K respectively for B = 1.0 X 10 12 , 3.0 x 10 12 and 
10 13 G. These have negligible effect on the acceleration but 
determine the initial ion charge Z. Actual values of Z a are 
unknown ah initio but are given by Z a — Z/(l+K), where K 
is the time-average of K at a specific position u on the polar 
cap and Z is here the pre-shower nuclear charge. Hence we 
assume arbitrary values Z a = 10 or 20 with mass numbers 
A = 20 or 40 and note that ions with very small Z s tend 
to be completely ionized by the polar-cap or whole-surface 
temperature at the start of acceleration and so produce no 
reverse-electron flux. The initial ion Lorentz factor is fixed 
as 7h = 5 at an altitude z — h = 0.05-R. The initial ion 
charges at z = h are: for Z s = 10; Zh = 8, 6 or 5 and for 
Z s = 20; Z h = 15, 12 or 9 respectively at B = 1.0 x 10 12 , 
3.0 x 10 12 G or 10 13 G.. 

Values of e a and Zoo obtained on this basis are given in 
Table 1. With the approximations and methods of Paper III 
in mind, it is obvious that the absolute values of e a do not 
have the accuracy that the number of figures given might 
indicate. But comparison of columns for different values of 
T a shows that whole-surface photons can generate very high 
transition rates not only at large V c for T a = 10 5 K, but 
for progressively lower V c as whole-surface temperature in- 
creases. By combining the results of Table 1, which are esti- 
mates of e a (V) and ZooiV), with equation (19) it should be 
possible, given Zh, to solve for V(u, oo). It can be seen by 
inspection of the equation, with reference to values of K (V) 
found from the Table, that there is always a solution, the 
value of V depending principally on T a and to a lesser extent 
on B. Values V <C Vmax occur at high T s and small B. 

Basically, incomplete ionization at z = h means that 
there will be further ionization at z > h unless T„ or Vmax, 
or both, are too small for it to be possible. In this latter 
case, relevant for old pulsars, there remains the polar-cap 
contribution e^ giving a small factor K < 1 and a solu- 
tion V(u, z) — V m ax(ti, z). Such a solution would be time- 
independent and of uniform p over the polar cap, with both 
ion and proton components. But it would not lead to ra- 
dio emission unless Vmax were small enough to give an ad- 
equate growth rate for the Langmuir mode, as described in 
Section 3. The general cut-off BP" 2 = 2.2 x 10 11 G s~ 2 ob- 
tained from the ATNF pulsar catalogue gives Knm( u i °°) = 
275(1 — u 2 /uq) GeV. Ion Lorentz factors 7 ~ 20 and pro- 
ton energies ~ 60 GeV are possible on flux lines originating 
near the edge of the polar cap. However, in the more gen- 
eral case, large values of e 3 in the uniform time-independent 



model lead to K S> 1 and hence to instability, as shown 
in Section 4.3 of Paper I. The actual state of the polar-cap 
surface is therefore time-dependent. 



5 QUALITATIVE MODELLING OF THE 
POLAR CAP 

The basic principle is that the state of any element Su of 
surface at any instant is a function only of its past time- 
distribution of reverse-electron energy flux. The flux of ac- 
celerated protons J p (u, t) is given by, 

J p (u,t) + J p (u,t) = f f p (t-t')K(u,t')J Z (u,t'), (20) 

in terms of the ion flux J z . The quantity J p — within 
intervals for which J p does not exceed the Goldreich- Julian 
current density paj(0)c- Most protons are produced in the 
vicinity of the shower maximum whose depth is only a slowly 
varying function of the primary electron energy. Their mo- 
tion toward the surface is better described by a drift veloc- 
ity, rather than diffusion, in most of the LTE atmosphere 
so that the distribution of arrival times can be modelled 
as f p (t — t') — 8(t — t' — t p ). We refer to Section 5 of 
Paper II for further details. Following Paper III, equation 
(20) recognizes that K is a function of u and t owing to 
its V — dependence, and can vary on time-scales of the or- 
der of t p . The protons are preferentially accelerated, but if 
J p reaches the Goldreich- Julian flux, as can be the case if 
K 2> 1, a thin proton atmosphere forms at the top of the 
ion atmosphere at a rate given by J p and the local ion flux 
falls to zero until it is exhausted. 

As a first approximation, equation (20) shows that the 
condition of an element of polar-cap area can be represented 
by two possible states: ion emission with duration t p or pro- 
ton emission of duration r gap ~ Kt p , where K is here the 
time-average over the interval t p . Values K 3> 1 ensure that 
proton production is so large that an atmosphere of protons 
forms at the top of the LTE atmosphere of ions. Then the 
ion phase ends and the proton phase is maintained until the 
proton component of the atmosphere is exhausted. Areas 
in the ion phase necessarily move with time and so can be 
seen as moving subpulses within the time-averaged radio- 
pulse profile. The motion could be chaotic or, as proposed 
in Paper II, organized so as to give regular subpulse drift. 

Pair creation by curvature radiation (CR) photons was 
referred to briefly at the end of Section 2 but has been oth- 
erwise ignored here. The necessary conditions are thought 
to be present in a relatively small fraction of the observed 
population (Harding & Muslimov 2002) but we shall review 
briefly its effects in the ft ■ B < case with SCLF bound- 
ary conditions. Self-sustaining pair creation requires ade- 
quate conversion probabilities for both outward and inward- 
directed CR photons at z < z max . Thus pair creation is not 
uniform over the whole cross-sectional area iru (z) of the 
open flux tube, but depends on a combination of V and of 
flux-line curvature. For the suitable values of u, the equa- 
tion analogous to (19), obtained by comparison of current 
densities at z = and 2 = z max , can be most compactly 
expressed in terms of the particle flux <j> on the relevant flux 
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lines, 

<t>P + l<?M = 1 | K f V(u,Oo) \ ^ 

4>p — \4>e\ 1-K \ Knai(u, OO) y ' 

in which <^> e is the flux either of electrons or positrons, and 
<f)p is the proton flux. This assumes, as is almost certainly 
the case, that the reverse-electron energy flux produces an 
excess of protons at the top of the LTE atmosphere. There 
is then no ion flow except on flux lines that do not satisfy 
the CR criteria such as those near ito on which Vmax is too 
small. 

Quantitative model construction is not attempted here. 
Owing to the uncertainties in significant parameters that are 
listed in Section 6, quantitative testing of a model against 
observed data does not appear immediately profitable. The 
purpose of this paper is less ambitious. It is an attempt to 
show that the condition of the polar-cap surface provides a 
physical basis for nulls, mode-switches and subpulse drift. 
To do this, we describe qualitatively some of the possible 
polar-cap states. 

Assume, initially, that CR pair production as described 
above is possible over some limited area of the polar cap. 
Its boundary can be defined as that within which proton 
production by the reverse-electrons of CR pairs exceeds the 
Goldreich- Julian rate so that a proton atmosphere forms 
and grows in density. Between this boundary and uo there 
is an annular strip within which ion acceleration can occur 
and proton formation is governed by equation (20). (Any 
contribution from CR reverse-electrons is neglected here.) 
Solutions of equation (20) then exist in which zones of ion 
phase circulate around the CR region. Equation (20) is local 
in u and so cannot itself determine the circulation sense or 
velocity. Suppose that there are n ion zones. Then in the 
usual notation by which subpulse drift is described, the cir- 
culation time is P3 = nPz, where P3 is the band separation 
and in the model is given by P3 = T gap + T v — (K + 1)t p . As- 
sociation with subpulse drift would require either secondary 
pair creation on the ion phase flux lines or the presence 
of proton and ion components in the current density near 
to uo so that energies are low enough for Langmuir-mode 
growth. A further consideration is that excess proton pro- 
duction within the CR area must be accompanied by lateral 
diffusion on the polar-cap surface, perpendicular to B. We 
have chosen the example of a clearly defined annular region 
but there is no reason why organized motion of this kind 
should not occur more generally and in the absence of a CR 
component. 

A further possibility is that the active CR region could 
be switched off by changes in the potential inside it which 
is also dependent on the state of the remaining parts of 
the polar cap. In general, in a less-organized state, there is 
no reason why either the whole polar cap or large parts of 
it should not be in the proton phase for intervals of time 
of the order of T gap and, in the absence of CR pair pro- 
duction, there would be no possibility of Langmuir mode 
growth. Thus nulls are a natural phenomenon within the 
model. Obviously, there are many possible polar-cap states, 
their relevance to widely-observed phenomena being limited 
only by the requirement that they should not be too depen- 
dent on neutron-star parameters having critical or particular 
values. 

Mode-changes have been studied extensively in a small 



Table 1. Values of the mean reverse-electron energy per ion (e s ) 
at the polar-cap surface in units of GeV, and of the mean ion 
charge at the end of photo-ionization (Zoo) have been calculated. 
The magnetic flux density is in units of 10 12 G and the accel- 
eration potential cut-off V c is in units of GeV. Arbitrary values 
Z a = 10 or 20 have been assumed for the mean surface nuclear 
charge. We refer to Section 4 for discussion of Z s and other pa- 
rameters on which the calculated values are dependent. The final 
three pairs of columns give the values of e s and of Z x for whole 
surface temperatures T s = 1, 2 and 4x 10 5 K. The rotation period 
is P = 1 s. Blank spaces denote energies that are negligibly small. 
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number of pulsars, some of which also exhibit subpulse drift. 
The paper of Bartel et al (1982) summarizes the main char- 
acteristics of this phenomenon. These authors emphasize 
that the spectral and polarization changes that occur on 
mode-switching are consistent with a change in the optical 
depth of the source region within the magnetosphere and 
that the sources of the two or more modes would be on differ- 
ent sets of flux lines. (They also suggest that the explanation 
of mode-switching lies in the polar-cap surface, though in the 
context of the Ruderman & Sutherland model.) Switching 
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of emission between different sets of flux lines occurs nat- 
urally in our qualitative model which is sufficiently diverse 
to accommodate the characteristics of mode-switching that 
are observed, including changes in subpulse drift rates. 

The qualitative model is applicable to subpulse drift 
and to short-duration nulls or mode changes and there 
is no difficulty in understanding how non-randomness or 
quasi-periodicity, such as that reported by Rankin & Wright 
(2008) and Redman & Rankin (2009) might occur naturally 
in such a system. The time-scales are necessarily short, per- 
haps of the order of T gap , but there is also the possibility 
of a medium time-scale instability in the atomic number Z s 
of surface nuclei, described in Section 4 of Paper II. This 
can result in values of Z 3 that are too small to produce 
any reverse-electron flux because the atoms are completely 
ionized in the LTE atmosphere. Having the largest charge- 
to-mass ratio, apart from protons, they are in equilibrium 
at the top of the atmosphere and so are preferentially accel- 
erated through the full potential difference V(u, oo) with no 
possibility of Langmuir-mode growth if CR pair creation is 
not possible. A broad measure of the time-scale for this in- 
stability is provided by the time in which a depth of surface 
nuclei equivalent to one radiation length is removed from 
the polar cap at the Goldreich- Julian current density. From 
equations (4) and (5) of Paper I, this is given by, 

tri=2 .lxl0 8 f " PS6C ^ 1/2 ) s. (22) 

V^Bi 2 ln(12Z 1 /2 J B 1 - 1 / 2 ); 

The distribution of Z 3 over the polar cap at any instant need 
not, and probably will not, be uniform. Thus complex dis- 
tributions of potential and of accelerated beam composition 
are to be expected. 



6 CONCLUSIONS AND UNCERTAINTIES 

Both neutron-star spin directions, defined by the sign of f2B 
at the polar caps, presumably exist and it is one of the pur- 
poses of this paper to question how, if at all, each sign con- 
tributes to the observed radio pulsar population. Previous 
papers (I - III) have been directed toward the CI • B < case 
because it has positive polar-cap corotational charge density 
so that the flux of accelerated particles can have ionic, pro- 
ton and positron components. There appears to have been 
some reluctance in the published literature to consider in 
any detail the physics of the polar cap for this spin direc- 
tion. But Papers I - III have attempted to show that such 
considerations are essential and that the radio emission char- 
acteristics are not solely determined by electrodynamics, as 
may be so for ft ■ B > 0. 

Weltevrede, Edwards & Stappers (2006) detected sub- 
pulse drift in almost one half of a sample of 187 pulsars se- 
lected only by signal-to-noise ratio. Nulls and mode changes 
are less frequently observed. We propose that results ob- 
tained in Papers I - III and in the present paper lead to 
the idea that these phenomena are all related to the physics 
of the polar cap in Q ■ B < pulsars and that the oppo- 
site spin direction, with only electrons as primary acceler- 
ated particles, is then observed as the population showing 
subpulse modulation, perhaps as a consequence of plasma 
turbulence, but none of the above phenomena. The problem 



with polar-cap physics is that a real condensed-matter sys- 
tem can be extremely complex and, as we noted in Section 
1, does not always lend itself to the construction of simple 
models (see also Section 2 of Paper II). A further difficulty is 
that some of the parameters concerned are either unknown, 
such as the atomic number of surface nuclei, or can vary 
over several orders of magnitude. 

Calculated values in Table 1 are limited to B = 10 12 — 
10 13 G which is representative of pulsars in the nulls tab- 
ulation of Wang et al (2007) and, for example, of those 
pulsars with drifting subpulses investigated by Gil et al 
(2008). Photo-electric cross-sections obtained in Paper III 
are not valid at B < 10 12 G in which region interpolation 
between zero-field cross-sections and those at 10 12 G would 
be necessary. At higher fields, B > B c = 4.41 x 10 13 G, 
approximate cross-sections for bremsstrahlung and pair cre- 
ation have been found (Jones 2010b) and electromagnetic 
shower development investigated. It was confirmed that the 
value of the proton production parameter W p decreases only 
slowly at B > B c but also that the Landau-Pomeranchuk- 
Migdal effect, known to exist in zero field, is significant in 
the neutron-star atmosphere, and has the effect of reducing 
high-energy bremsstrahlung and pair creation cross-sections 
and therefore of increasing shower depth and hence the pro- 
ton diffusion-drift time r p . 

The whole-surface temperature T s is a further impor- 
tant parameter that is not well known. Paper III did not in- 
corporate general-relativistic corrections that should strictly 
have been made in transforming the blackbody radiation 
field to the ion rest frame. But for most transitions, the ra- 
diation field is better approximated by that of the neutron- 
star surface proper frame rather than by T s °°, the temper- 
ature inferred by a distant observer. Since T s °° « 0.8T S for 
the neutron-star mass and radius assumed here in Section 2, 
the temperatures used in Table 1 are well below those that 
can presently be observed. Yakovlev & Pethick (2004) list 
only a small number of young pulsars whose blackbody tem- 
peratures have been measured. These are T%° > 5 x 10 5 K. 
All model predictions reviewed by these authors show T s °° 
decreasing very rapidly at age > 10 6 yr. But on the other 
hand, temperatures below T s = 10 5 K, the lowest used in 
Table 1, could be maintained by a very small and obscure 
level of dissipation within the star. Isolated neutron stars 
(INS) certainly have higher temperatures T s °° ~ 10 B K (see, 
for example, Kaplan & van Kerkwijk 2009). But they are 
few in number and their radio beams (if they exist) would 
be narrow compared with the 4-zr observability of unpulsed 
blackbody radiation. Given this uncertainty, we do not re- 
gard them as part of the population considered here. 

The energy of the coherent radio emission must pre- 
sumably be transferred from particle beams in high mag- 
netic field regions near the polar caps. Whatever the co- 
herent emission mechanism is eventually established to be, 
provided it is not curvature radiation, the criteria in Sec- 
tion 3 for growth of Langmuir modes must be relevant. This 
paper does not attempt to judge whether low-energy sec- 
ondary electrons or small Lorentz factor ions form the parti- 
cle beams producing radio emission in any particular group, 
but asserts that both possibilities exist for fi-B < pulsars. 
It follows that evolution of a pulsar to old age may include 
long intervals, perhaps of the order of 10 6 yr, in which pair 
creation is not possible because Vmax is too small but ion 
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Lorentz factors remain too large to give adequate Langmuir- 
mode growth rates. Because V max (u, oo) oc u 2 — it 2 , at least 
approximately, the highest growth rates will be on flux lines 
with u near uo- This is entirely consistent with the conal 
pattern of emission in the pulsar morphology developed by 
Rankin (1983) . 

Pulsars exhibit nulls with mean duration as short as 
several rotation periods or, for those described as intermit- 
tent pulsars, times of the order of 10 d. PSR 1931+24 is an 
example of the latter (Kramer et al 2006), but long-duration 
nulls have now also been observed in J 1832+0029 (see Lyne 
2009) and in J1841-0500 (Camilo et al 2012). Such null du- 
rations allow separate measurement of the spin-down rate in 
both on and off-states of emission. In each case, the off-state 
spin-down rate is about half that of the on-state, indicating 
some substantial difference in the condition of the magne- 
tosphere. Obviously, it is not known if this is also true for 
short-duration nulls but for fi ■ B < pulsars in our model of 
the polar cap, it occurs naturally owing to the change in the 
extent to which the electromagnetic fields near and beyond 
the light cylinder are loaded and the total momentum den- 
sity modified by particle emission. If the whole polar cap 
were in the proton phase, as discussed in Section 5, there 
would be no significant pair creation (assuming no CR pair 
production) and, of course, no possibility of Langmuir-mode 
growth. However, if creation of low-energy secondary pairs 
were possible, the flux of particles to be accelerated at or be- 
yond the light cylinder could greatly exceed the Goldreich- 
Julian value. Even though we lack a comprehensive under- 
standing of acceleration at Rlc and beyond, it is entirely 
plausible that the changes in particle number density in- 
herent in our model would lead to the observed changes in 
spin-down torque. 

Interesting reviews of pulsars with large nulling frac- 
tions and of the Rotating Radio Transients (RRAT) have 
been given recently by Burke-Spolaor & Bailes (2010), 
Keane et al (2011) and by Keane & McLaughlin (2011). 
In particular, pulsars are shown in the g — P plane, where g 
is here the fraction of periods in which a pulse is detected. 
The distribution appears roughly uniform in log g within the 
interval 10~ 4 < g < 10 _1 . The existence of these objects is 
not inconsistent with the complexity of possible polar-cap 
states of the fl-B < pulsars described here given the time- 
scale t r i defined by equation (22) . But Keane & McLaughlin 
note that there is a scarcity of objects whose behaviour is 
governed by time-scales that are longer than the above but 
shorter than the very long time-scales (~ 10 6 — 10 7 s) asso- 
ciated with the intermittent pulsars such as PSR 1931+24. 
Very long time-scales much exceeding those given by equa- 
tion (22) still remain a problem whose explanation may be 
more subtle than those attempted in this paper or may lie 
in quite different considerations. 

A simple explanation for the RRAT has been proposed 
by Weltevrede et al (2006). These authors compared RRAT 
pulses with the infrequent very large-amplitude pulses seen 
in PSR 0656+14. They note that if 0656+14 were as dis- 
tant in the galaxy as most of the RRAT, its observed emis- 
sion would appear to be that of the typical RRAT. But 
this may be inconsistent with the periodicities in individual 
RRAT pulse arrival times found recently by Palliyaguru et al 
(2011). These are mostly of the order of 10 4 s but a minority 
are much longer, of the order of 10 6 — 10 r s. If these com- 



plex sets of periodicities are established and are not loose 
quasi-periodicities, they represent a very real problem for 
the model of the polar cap described here. It is difficult to 
see how they can be present in a truly isolated neutron star, 
with either spin direction, not interacting with any other 
periodic system. Quasi-periodicities might be understand- 
able in terms of natural phenomena at the condensed- matter 
surface, possibly connected with the diffusion of \ow-Z a nu- 
clei perpendicular to B in the vicinity of the polar cap or 
with the two sections of the open magnetosphere described 
in Section 2. But well-defined periodicities would remain a 
problem. 
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